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We discuss renormalization effects on neutrino masses and mixing angles in a supersymmetric string- 
inspired SU(4)xSU(2)_LxSU(2)_RxU(f)x model, with matter in fundamental and antisymmetric 
tensor representations and singlet Higgs fields charged under the anomalous U(l)x family symmetry. 
The quark, lepton and neutrino Yukawa matrices are distinguished by different Clebsch-Gordan 
coefficients. The presence of a second U(l)x breaking singlet with fractional charge allows a more 
realistic, hierarchical light neutrino mass spectrum with bi-large mixing. By numerical investigation 
we find a region in the model parameter space where the neutrino mass-squared differences and 
mixing angles at low energy are consistent with experimental data. 



The experimentally measured values of gauge coupling 
constants a?, , a em , and the weak mixing angle sin 2 9w are 
correctly predicted in the Minimal Supersymmetric Stan- 
dard Model (MSSM) , assuming a unification scale of the 
order 10 16 GeV. Moreover, the existing data from neu- 
trino oscillation experiments provide an important clue 
to physics beyond the successful Standard Model (SM) 
and MSSM. 

Experimental data on atmospheric and solar neutrino 
oscillations 1] imply tiny but non-zero neutrino mass- 
squared differences Am 2 ... The negligible size of the 
neutrino masses, as compared to those of quarks and 
charged leptons, might suggest that a theory beyond 
MSSM should incorporate the right-handed neutrinos 
and an appropriate (see-saw) mechanism to suppress ad- 
equately the neutrino masses. Moreover, the observed 
large neutrino mixing angles 9 Ul2 , #„ 23 present challenges 
for additional symmetries and a unified framework in 
which neutrinos and quarks form part of same multiplet. 
Examples of higher symmetries including the SM gauge 
group and incorporating the right-handed neutrino in the 
fcrmion spectrum, are the partially unified Pati-Salam 
model, based on 577(4) x SU(2) L x SU(2) R 0, and the 
fully unified SO(10). When embedded into perturba- 
tive string or D-brane models, these may be extended by 
additional abelian or discrete fermion family symmetries. 
Thus fermion masses and mixing angles can be compared 
to the predictions of various types of models with full or 
partial gauge unification and flavor symmetries. 

Recently some models have been proposed to explain 
the presence of large mixing angles in the neutrino sec- 
tor, in contrast to the smaller quark mixings. For exam- 
ple, the mixing angle 9 Vl2 and the Cabibbo mixing 9c 
could satisfy the so called Quark-Lepton Complemen- 
tarity (QLC) relation 6 Ul2 + 8 C ~ f Q- It has been 
shown Q that this relation can be reproduced if some 
symmetry exists among quarks and leptons. Attempts 
to realize QLC in the context of models unifying quarks 
and leptons such as Pati-Salam have been made y|. 

As another possibility, the symmetry L e — L^ — L T im- 
plies an inverted neutrino mass hierarchy and bimaximal 



mixing 9 Vl2 = 9„ 23 = f, with 0„ 13 =0 6]. This symme- 
try alone does not give a consistent description of current 
experimental data, but additional corrections and renor- 
malization effects have still to be taken into account. It 
has been shown Q in the context of MSSM extended by 
a spontaneously broken U(l)x factor, that the neutrino 
sector respects an L e — — L T symmetry. Small correc- 
tions from other singlet vevs, which are usually present 
in a string spectrum, can lead to a soft breaking of this 
symmetry and describe accurately the experimental neu- 
trino data. 

Another important issue is the renormalization group 
(RG) flow of the neutrino parameters from the high en- 
ergy scale where the neutrino mass matrices are formed, 
down to their low energy measured values. One can at- 
tempt to determine the neutrino mass matrices from ex- 
perimental data directly at the weak scale. However, the 
Yukawa couplings and other relevant parameters are not 
known at the unification scale. A knowledge of these 
quantities at the unification mass could provide a clue 
for the structure of the unified or partially unified the- 
ory and the exact (family) symmetries determining the 
neutrino mass matrices at the GUT scale. Attempts to 
determine the neutrino mass parameters in a top-bottom 
or bottom-up approach have recently discussed in the lit- 
erature [III El. 

In this paper we investigate further the neutrino mass 
spectrum of a model with gauge symmetry SU(4) x 
SU(2) x SU{2) x U(l) x based on the 4-4-2 models 
Il4l | , whose implications for quark and lepton masses were 
recently investigated in ^lj. These models present sev- 
eral attractive features. Firstly, only "small" Higgs repre- 
sentations are needed and these commonly arise in string 
models. Secondly, third generation fermion Yukawa cou- 
plings are unified [l5| up to small corrections. Unifica- 
tion of gauge couplings is allowed and, if one assumes 
the model embedded in supersymmetric string, might be 
predicted Furthermore, the doublet-triplet splitting 
problem is absent. 

In string derived models a large number of neutral sin- 
glet fields carrying quantum numbers only under U(l)x 
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appear in the spectrum of the effective field theory model. 
D- and F-flatness conditions require some of them to ob- 
tain non-zero vevs of the order of the \J(X)x breaking 
scale. In the present model, in order to describe accu- 
rately the low energy neutrino data we introduce two 
such singlets charged under U(l)x [lq - The previous 
model [TlJ with one such singlet could easily give rise to 
a spectrum of light neutrinos with normal hierarchy and 
bi-large mixing. However, after study of the renormaliza- 
tion group (RG) evolution and unification it was found 
that the scale of light neutrino masses too large to be 
compatible with observation. If we impose the correct 
scale of light neutrino masses, then some heavy right- 
handed neutrinos (RHN) would have masses above the 
unification scale, which is incompatible with our effec- 
tive field theory approach. 

Thus three a priori independent expansion parame- 
ters arise from the superpotential, two coming from the 
singlets and one from the Higgses which receive v.e.v.'s 
at the SU(4)xSU(2)/{ breaking scale. In general, a 
nonrenormalizable operator may contain several prod- 
ucts of the SU(4)xSU(2)fl breaking Higgses, and thus 
different contractions of gauge group indices are possi- 
ble leading to different contributions depending on the 
Clebsch factors. We use a minimal set of nontrivial 
Clebsch factors to construct the Dirac mass matrices. 
Right-handed (SU(2)l singlet) neutrinos acquire Majo- 
rana masses through nonrenormalizable couplings to the 
U(l)x-charged singlets and to Higgses, while light neu- 
trinos will obtain masses via the see-saw mechanism. 

The renormalization group equations (RGEs) for the 
neutrino masses and mixing angles above, between and 
below the see-saw scales are solved numerically, for sev- 
eral sets of order 1 parameters which specify the heavy 
RHN matrix. In each case the results at low energy are 
consistent with current experimental data, and provide 
further predictions for the 1-3 neutrino mixing angle and 
for neutrinoless double beta decay. 



which acquire vevs of the order Mq along their neutral 
components 



(H) = (v H )=M G , (H) = {u R ) = M c 



(1) 



The Higgs sector also includes the h = (1, 2, 2) field which 
after the breaking of the PS symmetry is decomposed 
to the two Higgs superfields of the MSSM. Further, two 
D = (6, 1, 1) scalar fields are introduced to give mass to 
color triplet components of H and H via the terms HHD 
and HHD O. 

Finally, we introduce two scalar singlet fields <fi, 
X, charged under U(l)x whose vevs will play a cru- 
cial role in the fermion mass matrices through non- 
renormalizable terms of the superpotential. In the sta- 
ble SUSY vacuum the two singlets obtain vevs to satisfy 
the D-flatness condition including the anomalous Fayet- 
Iliopoulos term [24[. The anomalous D-flatness condi- 
tions allow solutions where the vevs of the conjugate 
fields (j) an d X are zero and we will restrict our analy- 
sis to this case. Note that in general a string model may 
have more than two singlets and more than one set of 
Higgses H i: Hi, with different V(l)x charge. All such 
fields may in principle also obtain vevs, however we find 
that two of them are sufficient to give a set of mass ma- 
trices in accordance with all experimental data. Hence 
we consider any additional singlet vevs to be significantly 
smaller. 

The Higgses Hi, Hi may obtain masses through HHtfi, 
HHx and HH(f>x couplings. However, in order to break 
the Pati-Salam group while preserving SUSY we require 
that one H-H pair be massless at this level. This 
"symmetry-breaking" Higgs pair could be a linear combi- 
nation of fields with different U(l)x charges, which would 
in general complicate the expressions for fermion masses. 
The chiral spectrum is summarized in TableQ] We choose 
the charge of the Higgs field h to be —013 — so that 
that the 3rd. generation coupling F3F3I1 is allowed at 
tree-level. 



I. DESCRIPTION OF THE MODEL 

In this section we present salient features of the 
string inspired Pati-Salam model extended by a 
\J(l)x family-symmetry, the total gauge group being 
SU(4)xSU(2) L xSU(2) fl xU(l) x . The field content in- 
cludes three copies of (4,2, 1) + (4, 1,2) representations 
to accommodate the three fermion generations Fi + F L 
(4 = 1,2,3), 



Fi = 



Ui Ui 

di di 



Fi = 



Ui u\ u\ 
di df di 



where the subscripts on, cti indicate the U(l)jf charge. 
In order to break the Pati-Salam symmetry down to SM 
gauge group, Higgs fields H = (4, 1, 2) and H = (4, 1, 2) 
are introduced 
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TABLE I: Field content and U(l)x charges 
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We now turn to the terms in the superpotential which 
can give rise to fermion masses. Dirac type mass terms 
arise after electroweak symmetry-breaking from cou- 
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plings of the form 



II. FERMION MASS MATRICES 
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Apart from the heaviest generation, all masses arise at 
non-renormalizable level, suppressed by powers of the 
fundamental scale or unification scale Mtj. The couplings 
Vmm'' v'n e ^ c - are non- vanishing and generically of or- 
der 1 whenever the \J(l)x charge of the corresponding 
operator vanishes, thus: 

a.i — a 3 + otj — &3 — {—mz, —m'z', —n(x + x), 

— kz — £z , —p(x + x) — qz, —r(x + x) — sz'}. 

Other higher-dimension operators may arise by multiply- 
ing any term by factors such as (H H) e (j> s /M^ +s where 
l(x + x) + s z = 0. Such terms are negligible unless the 
leading term vanishes. 

Neutrinos may in addition receive also Majorana type 
masses. These arise from the operators 
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A. General structure 

As can be seen from the superpotential Yukawa cou- 
plings J2J and (0), three different expansion parameters 
appear in the construction of the fermion mass matrices. 
These are 



_ m 



Ml 



Ml 



Ml 



Mtj 



given 



(HH) 



M%. 



Note that, for non-renormalizable 



Dirac mass terms involving several products of HH/My, 
the gauge group indices may be contracted in different 
ways 141. This can lead to different contributions to 
the up, down quarks and charged leptons, depending on 
the Clebsch factors C^ u dev \ multiplying the effective 
Yukawa couplings. Also, although the Clebsch coeffi- 
cient for a particular operator O n may vanish at order 
n, the coefficient for the operator 0^ n+p y q containing p 
additional factors (HH) and q factors of 4> and/or \ is 
generically nonzero. 

In our analysis we wish to estimate the effects of the 
second singlet (x) contributions on the neutrino sector 
as compared to the analysis presented in |ll| without 
affecting essentially the results in the quark sector. In 
order to obtain a set of fermion mass matrices with the 
minimum number of new operators, we assume fractional 
U(l)x charges for H, H and x fields, while the combina- 
tion HH and the singlet <j> are assumed to have integer 
charges. Thus a*, oti, x + x and z are integers, while z', 
x and x are fractional. As a result, the Dirac mass terms 
involving vevs of x are expected to be subleading com- 
pared to other terms. Suppressing higher-order terms 
involving products of e, e' and e", the Dirac mass terms 
at the unification scale are 

rtH, « S a 6 j3 m 3 + (e m + (e") m ' + C y (e')") v u , d (5) 

where 7713 = v n .<iyo 3 , with v u and Vd being the up- type 
and down-type Higgs vevs respectively, and we omit the 
order-one Yukawa coefficients etc. for simplicity. 

The Majorana mass terms are proportional to the com- 
bination HH (see Eq. ©) which has fractional U(\)x 
charge. Thus, terms proportional to x/Mjj become now 
important for the structure of the mass matrix. The gen- 
eral form of the Majorana mass matrix is then 



Couplings of this type are non-vanishing whenever the 
following conditions are satisfied: 



on + aj + 2x — {— tz, —tz, —w(x + x), 

- k'z - t'z' - p'{x + x) - q'z, -r'(x + x) - s'z'}. 



M N « Mr (a4V + ${e"f + 



where we define Mr = Mq/Mtj = e'Mrj. 
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B. Specific choice of U(l)x charges 

Before we proceed to a specific, viable set of mass ma- 
trices, we first make use of the observation that the 
form of the fermion mass terms above is invariant under 
the shifts 
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TABLE II: Possible choices for the U(l)x charges of the H 
and x fields. 



+ C, a-i — > a« + C £ 



C, a; 



C (6) 



so that we are free to assign as, = az = 0. We further fix 
x + x — 1 and z = — 1; we will choose the values of x and 
z' to be fractional such that the v.e.v. of \ only affects 
the overall scale of neutrino masses, as explained below. 

The charge entries of the common Dirac mass matrix 
for quarks, charged leptons and neutrinos are then 



Qx[M. 



D\ 




(7) 



and the charge matrix for heavy neutrino Majorana 
masses is 



Qx[M N ] = 2x 




(8) 



Now, we relate e, e', e" with a single expansion parameter 
7?, assuming the relations 



Vv 



(9) 



where b\ , 62 are numerical coefficients of order one. Then 
the effective Yukawa couplings for quarks and leptons 
may include terms 

Yf = b lV m/2 + &™+y+™/2 + cjrf 1 ' 2 + b 1 r 1 1+n ^ 2 + ■■■ 

(10) 

with / = u,d, e, v, up to order 1 coefficients yV . Which of 
these terms survives, depends on the sign of the charge of 
the corresponding operator. For a negative charge entry, 
the first two terms are not allowed and only the third and 
fourth contribute. Further, if a particular CV coefficient 
is zero, then we consider only the fourth term. 

Therefore, we need to specify the Clebsch-Gordan coef- 
ficients CV for the terms involving powers of (HH) /My. 
These coefficients could be found if the fundamental the- 
ory was completely specified at the unification or string 
scale. In the absence of a specific string model, here we 
present a minimal number of operators which lead to a 
simple and viable set of mass matrices. Up to possible 
complex phases, we choose C\? = Cj 2 = |, C 23 = 3 and 
Cl 1 = = Cf = Cf = Cf = C 22 = Cf = with 
all others being equal to unity. The effective Yukawa 
matrices at the GUT scale obtained under the above as- 



sumptions are 
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(11) 



where we suppress order one coefficients. The quark sec- 
tor as well as the neutrino sector were studied in 
However, full renormalization group effects were not cal- 
culated for the neutrino sector and as it turns out one 
singlet is inadequate to accommodate the low energy 
data. Consequently, we introduced the second singlet 
X, with fractional charge, whose v.e.v. affects only the 
overall scale of neutrino masses. 

The desired matrix for the right handed Majorana neu- 
trinos may result from more than one choice of charge for 
the H field and the x singlet field. These can be seen in 
Table HT1 We choose the H charge to be x : 
2 x is non-integer, and set the x singlet charge to 
The analysis for the quarks and charged leptons remains 
the as in since operators with nonzero powers x r do 
not exist for powers r < 5 and are negligible compared 
to the leading terms. 

With these assignments, the charge entries of the heavy 
Majorana matrix Eq. © are: 
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(12) 



Due to the fractional U{X)x charges contributions from 
(j> or HH alone vanish. However, we also have the sin- 
glets HHx/Mjj with vev 6 2 rf^ 2 an, l 4>x/Mu a vev 
^1^2 f? 2 , while for some entries one may have to consider 
higher order terms since the leading order will be van- 
ishing. In Table IIIII we explicitly write the operator for 
every entry of M/v- The Majorana right-handed neutrino 
mass matrix is then 
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(13) 



with M R = e'Mu = y/rjMu. 

Having defined the Dirac and heavy Majorana mass 
matrices for the neutrinos, it is straightforward to obtain 
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TABLE III: Operators producing the Majorana right handed 
neutrino matrix Mn- 

the light Majorana mass matrix from the see-saw formula 

m„ = -m Dl/ M N 1 m T DlJ (14) 
at the GUT scale. 

C. Setting the expansion parameters 

Given the fermion mass textures in terms of the U (l)x 
charges and expansion parameters, we need now to de- 
termine the values of the latter in order to obtain consis- 
tency with the low energy experimentally known quanti- 
ties (masses and mixing angles). Note that the coefficient 
62 defined in Eq. © will determine the overall neutrino 
mass scale through Eq. (|13[1 . 

Consistency with the measured values of quark masses 
and mixings fixes the value of 77 « 5 x 1CP 2 : for example 
the CKM mixing angle #12 is given by y/rj ps 0.22 up 
to small corrections [ljj. Hence the ratio of the SU(4) 
breaking scale Mq to the fundamental scale Mjj is also 

A/ 2 

fixed through -^gf- = ^/rj w 0.22: the Pati-Salam group is 
unbroken over only a small range of energy. We perform 
a renormalization group analysis in order to check the 
consistency of this prediction with the low-energy values 
of the gauge couplings a s , a em and the weak mixing angle 
sin 2 9w [l3 

sin 2 9 W = 0.23120, a 3 = 0.118 ± 0.003, a em = 

127.906 

If the underlying model at Mjj has a single unified gauge 
coupling, then Mq is fixed to be just below the unifi- 
cation scale according to the analysis of gauge coupling 
unification in the MSSM. Because of this fact, the low 
energy measured range for as affects the unification of 
the gauge couplings. Thus, we add the following extra 
states 

h L = (1,2,1), h R = (1,1,2) 

which are usually present in a string spectrum |12| . It 
turns out that we need 4 of each of these extra states for 



Mq = 9.32 x 10 15 GeV to be consistent with the value 
of e' deduced from quark mass matrices. 

In Figure^we plot the evolution of the gauge couplings 
from Mz to Mjj. In Figure [21 we show in more detail the 
evolution of the gauge couplings in the Pati-Salam energy 
region. The two bands for the 0:4 and ctiR couplings are 
due to strong coupling uncertainty at Mz- For a%{Mz) = 
0.1176, as can be seen from Figure |3 we obtain Mjj = 
1.96 x 10 16 GeV. 
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FIG. 1: Evolution of the gauge couplings. The two lines for 
Q!3 indicate the range of initial conditions at Mz- 
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FIG. 2: Close-up of the gauge couplings in the Pati-Salam 
energy region. 

The remaining parameters to be determined are &i , 62 
and fiij of the right handed neutrino mass matrix. We 
find that Mn is proportional to &2, Mn « b2MnM' N , 
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FIG. 3: The unification point of Pati-Salam gauge couplings 
for a 3 (M z ) = 0.1176. 



thus 62 is related to the scale of the light matrix m„. 
Also, the choice b\ w 1.1 leads to agreement with the 
data, while implying e = ~ 0.25. 



III. RUNNING OF NEUTRINO MASSES AND 
MIXING ANGLES 

One of the problems one encounters when searching for 
a specific mass matrix for the light neutrinos via the see- 
saw formula is the effects induced by the renormalization 
group equations. The low energy neutrino data could 
be considerably different from the results at the see-saw 
scale. The running of neutrino masses and mixing angles 
has been extensivel y di scussed for energies below the see- 
saw scales [H as well as above H H H^. The 
running of the effective neutrino mass matrix m v above 
and between the see-saw scales is split into two terms, 

m v = -— ( K + 2Y l ,M- 1 Yj ) . (15) 

where K is related to the coefficient of the effective 5 di- 
mensional operator LLh u h u , (n) labels the effective field 
theories with M n right handed neutrino integrated out 

(M n > M n _i > M„_2, . . .) and Y v are the neutrino cou- 
plings at an energy scale M between two RH neutrino 

(«) 

masses M n > M > M n _i, while Y u = below the light- 
est RH neutrino mass. These effective parameters govern 
the evolution below the highest seesaw scale and obey the 



differential equation [TIIMI^ 

M 

167T 2 — = (F e r e t + Y u Yj) T X + X(Y e Y^ + Y u Yj) T 

+ (2Tt(yJyJ + ZY u Y^) - 6/5 3 2 - 6g 2 2 )X (16) 

where X = k, Y u M^ 1 Yj . The RGEs have been solved 
both numerically and also analytically H, 0, |2l| . Numer- 
ically, below the lightest heavy RH neutrino mass large 
renormalization effects can be experienced only in the 
case of degenerate light neutrino masses for very large 
tan f3 0, |2^| . Above this mass things are more com- 
plicated due to the non-trivial dependence of the heavy 
neutrino mass couplings, unless M v ? is diagonal. For the 
leptonic mixing angles, in the case of normal hierarchy 
relevant to our model, one expects negligible effects for 
the solar mixing angle while #13 and #23 are expected to 
run faster . 

On the other hand, studying the analytical expressions 
obtained after approximation, exactly the opposite be- 
havior is predicted and the solar mixing angle receives 
larger renormalization effects than #13 or #23- However, 
possible cancellations may occur and enhancement or 
suppresion factors may appear, thus the numerical so- 
lutions may differ considerably from these estimates. 

In our string-inspired model the Dirac and heavy 
Majorana mass matrices at the unification scale are 
parametrized in terms of order- 1 superpotential coeffi- 
cients fiij(Mu) whose exact numerical values are not 
known. The flavour structure at the unification scale 
might also be different from that at the electroweak scale 
Mz- Thus, even if the Yukawa parameters are deter- 
mined at Mz, to understand the structure of the theory 
at Mjj, and consequently any possible family symmetry, 
we would certainly need the parameter values at Mjj. 

In this section we study the renormalization group flow 
of the neutrino mass matrices "top-down" from the Pati- 
Salam scale Mq to the weak scale. We choose sets of val- 
ues of the undetermined order 1 coefficients at the high 
scale and run the renormalization group equations down 
to Mz where we calculate AmJ^ and 6 Uij and compare 
them with the experimental values. Study of a bottom- 
up approach has been performed ^} and we will compare 
our results to this work. The renormalization group anal- 
yses of the neutrino parameters, successively integrating 
out the right handed neutrinos, is performed using the 
software packages REAP/MPT @. 

• We generate appropriate numerical values for the 
coefficients /in, /112, /113, /122, M23, so that after the evo- 
lution of ra v to low energy we obtain values in agreement 
with the experimental data. The coefficient /i33 is set to 
unity (which can always be done by adjusting the value 
of b-i). Experimentally acceptable solutions can be seen 
in Table llVl In Table Ivl we present the resulting values of 
Oij and Ami at the scale Mz- The mass-squared differ- 
ences lie in the ranges Am^ fm = [1.33 — 3.39] x 10 _3 eV 2 , 
Am 2 sol = [7.24 - 8.85] x !CT 5 eV 2 . These are consistent 
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Solution 


/in 


/X12 


/Ul3 


(122 


M23 


1. 


0.10535 


0.10972 


0.86012 


0.10491 


0.91014 


2. 


0.11939 


0.10954 


0.80912 


0.10683 


0.93832 


3. 


0.10392 


0.11787 


0.97796 


0.10512 


0.98749 


4. 


0.09143 


0.10962 


0.87616 


0.10063 


0.93798 


5. 


0.12697 


0.12745 


0.99860 


0.11652 


0.99980 


6. 


0.10920 


0.09638 


1.00975 


0.10238 


0.93561 


7. 


0.10124 


0.11682 


0.98568 


0.10688 


0.99156 


8. 


0.12358 


0.09514 


0.99580 


0.10434 


0.95646 


9. 


0.13006 


0.11973 


1.02235 


0.10378 


0.89460 


10. 


0.12665 


0.12137 


1.00029 


0.10695 


0.91578 



> 

CD 




10 12 14 16 

log 10 (n/GeV) 



TABLE IV: Numerical values of parameters /Ltii, (112, M13, 
M22, M23 at Ma- 



FIG. 4: The running of the light neutrino masses. 



Soln. 


Am 2 atm ■ 10 3 


AmL ■ 10 5 


6*12 


013 


023 


1. 


2.7149 


7.9621 


29.442 


3.9859 


44.114 


2. 


2.3145 


7.9514 


34.289 


12.507 


51.047 


3. 


1.8978 


8.6141 


30.560 


0.8656 


46.230 


4. 


3.0062 


8.3217 


34.347 


1.8512 


44.333 


5. 


3.3905 


7.2468 


30.245 


2.9355 


36.900 


6. 


3.2459 


7.5351 


34.296 


1.3701 


46.947 


7. 


2.0171 


7.9464 


34.432 


1.0086 


50.279 


8. 


1.3321 


7.9060 


37.646 


6.1490 


43.067 


9. 


2.4867 


8.8561 


29.592 


5.6007 


42.970 


10. 


2.1652 


7.8869 


29.189 


3.1512 


37.220 



10 12 14 16 

log 10 (n/GeV) 



TABLE V: Values of the physical parameters Anv^ tm , Am 2 aol , 
0i2, 013, 023 at Mz (mass units eV 2 ). 



FIG. 5: The evolution of the mixing angles. 



with the experimental data Am^ 
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atm.exp ~ [1-3 3.4] X 

3 eV 2 and Am 2 solexp = [7.1 - 8.9] x lCT 5 eV 2 . The 
mixing angles are also found in the allowed ranges 6*12 = 
[29.4-37.6], 6» 23 = [36.9-51.0] and 13 = [0.86- 12.50]. 

• In figure (QJ we plot the running of the three light 
neutrino Majorana masses (mi < m 2 < to 3 ) in the en- 
ergy range Mq — Mz- The initial (GUT) neutrino eigen- 
masses are all larger than their low energy values. Signif- 
icant running is observed mainly for the heaviest eigen- 
mass . For experimentally acceptable mass-squared 
differences Am 2 ,. . at Mz, in all cases their corresponding 
values at the GUT scale lie out of the acceptable range. 
In this scenario with hierarchical light neutrino masses, 
we find that large renormalization effects occur above the 
heavy neutrino threshold since the Yukawa couplings Y v 
are large and the second term in (|15|) dominates. Also, 
since m Vl < y/ Am 2 o[ , the solar angle turns out to be 
more stable compared to the running of the 623, as can 
be seen in Figure 10. These results are in agreement 
with the findings of jicj . 

• In figure (0 we plot the distribution Am 2 tm versus 
Am 2 sol at the two scales M G (Table EJJ and M z for the 



ten models of Table ITVl We find that the hierarchy of the 
neutrino masses at the Pati-Salam breaking scale tends 
to be greater than that at low energies. Several mod- 
els predict Ato^/Ato^ out of the experimental range at 
Ma, although after the running at Mz they are consis- 
tent with the data. 

• Finally, we check the predictions of our model for 
the effective neutrino mass parameter relevant for Pf3o v 
decay. This parameter can be written in terms of the 
observable quantities as 



I (to) I = (mi cos 2 6> + e la V 'Am soi sin 2 9 Q ) cos 2 9 13 



TO Qtm sm «i 3 e 



(17) 



In the last column of Table l|VI|) the /3/3o^-decay pre- 
dictions are presented for solutions 1-10. Many current 
experiments attempt to measure this quantity [25J; the 
best current limit on the effective mass is given by the 
Heidelberg-Moscow collaboration [2(| 



(to) < 0.35 z eV, 



(18) 
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Soln. 


Am 2 atm (M G ) ■ 10 3 


AmL(Mc) ■ 10 5 


(m) 


hi 


1. 


10.5294 


27.6096 


0.00361 


3.41 


2. 


8.1987 


30.1333 


0.00633 


2.88 


3. 


7.2533 


31.6108 


0.00607 


1.50 


4. 


11.7749 


30.4352 


0.00753 


1.28 


5. 


14.093 


23.7416 


0.00346 


2.85 


6. 


12.347 


28.2642 


0.00822 


1.26 


7. 


7.4256 


30.855 


0.00821 


1.23 


8. 


5.2910 


29.8775 


0.00852 


1.17 


9. 


9.74656 


30.5418 


0.00379 


3.64 


10. 


8.99695 


25.8838 


0.00327 


3.42 



TABLE VI: Values of the physical parameters Am„, m and 
Ara^i at Mg\ the effective neutrino mass (m) related to PPo* 
decay; and the parameter t>2 which determines the scale of the 
light matrix m v . 



where the parameter z = 0(1) allows for uncertainty 
arising from nuclear matrix elements. 

In a recent analysis of neutrinoless double beta de- 
cay the allowable range of the effective mass param- 
eter was given for specific scenarios. In the case of the 
normal hierarchy the bounds are 



< (m) < 0.007 eV 



(19) 



thus our results are in the experimentally acceptable re- 
gion. 



E 

< 10"' 



10 12 14 16 

log 10 ( M /GeV) 



FIG. 6: Running of AmJ oJ and AmJ t 



™ -1.0 - 
E 
3* 

cn -1.1 - 
o 

-1.2 - 



At M ; 
At M r 



-2.05 -2.00 -1.95 -1.90 -1.85 -1.80 -1.75 

2,1/2 



log, (Am,;)" 
FIG. 7: 

^^12 cLnd A77123 at Mz and at Mg- 

constructed in terms of three expansion parameters 
_ (HH) 



M v ' 



Ml 



Mu 



(20) 



where 4> and \ are singlets and H , H the SU(4) x SU(2)fl- 
breaking Higgses. The model is simplified by the frac- 
tional U(l)x charges of H and \-> which ensure that the 
parameter e" only appears as a prefactor to the heavy 
Majorana neutrino masses. 

The expansion parameter arising from the Higgs 
v.e.v.'s defines the ratio of the £(7(4) breaking scale Mq 
to the unification scale Mjj: we performed a renormal- 
ization group analysis of gauge couplings under this con- 
straint and found successful unification with the addition 
of extra states usually present in a string spectrum. 

Assuming that only the third generation of quarks and 
charged leptons acquire masses at tree level and under a 
specific choice of U(l)x charges as well as Clebsch fac- 
tors, we examined the implications for the light neutrino 
masses resulting from the see-saw formula. We found 
that the light neutrino mass spectrum is hierarchical and 
that the mass hierarchy tends to be larger at the GUT 
scale than at Mz due the rcnormalization group running. 
The solar mixing angle 612 is stable under RG evolution 
while larger renormalization effects are found for the at- 
mospheric mixing angle 823 and #13, always with their 
values at Mz in agreement with experiment. 



IV. CONCLUSIONS 

In this work, we studied the running of neutrino masses 
and mixing angles in a supersymmetric string-inspired 
SU(4)xSU(2) L xSU(2) fl xU(l) x model. An accurate de- 
scription of the low energy neutrino data forced us to in- 
troduce two singlets charged under the U(l)x, leading to 
two expansion parameters. The mass matrices are then 
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